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ABSTRACT 

A new age-redshift test is proposed in order to constrain Ho with basis on 
the existence of old high redshift galaxies (OHRG). As should be expected, the 
estimates of H based on the OHRG are heavily dependent on the cosmological 
description. In the flat concordance model (ACDM), for example, the value of 
H depends on the mass density parameter Qm = 1 — Oa- Such a degeneracy 
can be broken trough a joint analysis involving the OHRG and baryon acoustic 
oscillation (BAO) signature. In the framework of the ACDM model our joint 
analysis yields a value of H = 71 ± 0.04 km s -1 Mpc -1 (la) with the best fit 
density parameter Q M = 0.27 ± 0.03. Such results are in good agreement with 
independent studies from the Hubble Space Telescope key project and the recent 
estimates of WMAP, thereby suggesting that the combination of these two in- 
dependent phenomena provides an interesting method to constrain the Hubble 
constant. 

Subject headings: Hubble constant, galaxy, age problem, baryon acoustic oscilla- 
tions 
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1. INTRODUCTION 

The determination of Hq has a practical and theoretical importance to many astro- 
physical properties of galaxies, clusters and quasars, and plays an important role for several 
cosmological calculations and others physical rulers, like SNe la, the present energy density 
scale, primordial nucleosynthesis, and the age of the Universe (Peebles 1993; Peacock 1999; 
Freedman 2000). However, Spergel et al. (2007) have shown that CMB studies can not 
supply strong constraints on the value of H Q on their own. This problem occurs due to 
the high degree of degeneracy on the parameter space (Tegmark et al 2004), and may be 
circumvented only by using independent measurements of H (Hu 2005, Cunha et al. 2007). 
In these circumstances, it is particularly relevant to obtain accurate and independent bounds 
on the value of Ho from physics relying on many different kinds of observations. 

On the other hand, with the new optical and infrared techniques together the recent 
advent of large telescopes became possible to estimate the ages of high redshift objects like 
galaxies and quasars. The subsequent discovery of some evolved high-z objects (Dunlop et 
al. 1996, Spinrad et al. 1997; Bruzual & Magris 1997; Dunlop 1999; Nolan et al. 2001) 
have been used by many authors to put limits on the cosmological parameters. Alcaniz and 
Lima (1999, 2001), by using only two old (extremely red) high redshift galaxies (53W091 
and 53W069) dated by Dunlop and collaborators obtained f^A > 0.5 for a flat ACDM model, 
a result more stringent than the ones derived from Type la Supernovae. They also derived 
limits on the cosmic equation of state parameter by considering a flat cosmology driven by 
cold dark matter and a X-matter dark energy component (p = up,u = constant). As a 
conclusion, it was argued that if the ages of those objects were correctly estimated, the first 
formation era is pushed back to extremely high redshifts. Hasinger, Schartel and Komossa 
(2002) also reported the discovery of the quasar APM 08729 + 5255 at z = 3.91 with an 
estimated age of 2-3 Gyr. Later on, this result was independently confirmed by Friaca et 
al. (2005) by using a chemodynamical model for the evolution of spheroids. They shown 
that a minimal age of 2.1 Gyr is set by the condition that the F e /0 abundance ratio of 
the model reaches 3.3 (normalized to solar values) which is the best fit values obtained by 
Hasinger and coworkers. Further, several authors have also claimed that the quasar APM 
08729 + 5255 is an extremely stringent object for any realistic cosmological model (Alcaniz, 
Lima & Cunha 2003; Cunha & Santos 2004, Deepak and Dev 2006). In this regard, it is of 
particular interest the sample of 20 old passive galaxies published by the team of the Gemini 
Deep Deep Survey (GDDS). As noticed by McCarthy et al. (2004), such data are strongly 
suggesting that the galaxies have been evolving passively since their initial burst of star 
formation. More recently, the GDDS sample has been used by Dantas et al. (2007) to put 
constraints on the dark energy equation of state (w) through a joint analysis involving the 
inferred lookback time measurements from the ages of the galaxies plus the total expanding 
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age of the Universe as obtained from CMB data. 

In this letter, we derive new constraints on the Hubble constant H , by using the GDDS 
sample complemented by two additional data set recently released (Roche et al. 2006, 
Longhetti et al. 2007). A statistical age test is performed to constrain the (h,Q m ) plane 
assuming that the Universe is well described by the cosmic concordance model (ACDM). Due 
to a degeneracy on the space parameter, more stringent constraints are obtained through a 
joint analysis involving the different OHRGs samples and the current SDSS measurements 
of the baryon acoustic peak (Eisenstein et al. 2005). In this concern, we recall that the 
baryon acoustic oscillations (BAO) method is independent of the Hubble constant H$ but 
it is heavily dependent on the value of Om, and, as such, it contributes indirectly to fix 
the value of H . Therefore, one can apply the BAO signature to break the degeneracy of 
the mass density parameter (f2 M ). The BAO signature comes out because the cosmo logical 
perturbations excite sound waves in the relativistic plasma, thereby producing the acoustic 
peaks in the early universe. Eisenstein et al. presented the large scale correlation function 
from the Sloan Digital Sky Survey (SDSS) showing clear evidence for the baryon acoustic 
peak at lOO/i -1 Mpc scale, which is in excellent agreement with the WMAP prediction from 
the CMB data. As we shall see, both set of data (OHRG sample and BAO signature) provide 
an interesting method to constrain the Hubble constant. 



2. Basic Equations and Sample 

Let us now consider that the Universe is described by a flat Friedmann-Robertson- 
Walker (FRW) geometry (c = 1) 

ds 2 = dt 2 - a 2 it) [dr 2 + r 2 (d9 2 + sin 2 0d<f>)] , (1) 

as theoretically predicted by inflation and required by the WMAP observations (Spergel et 
al. 2003; 2007). 

In such a background, the age-redshift relation for a model driven by cold dark matter 
plus a cosmological constant (ACDM) has only two free parameters (H ,Qm)- It can be 
written as 

r { : l+z ^ 1 dr 

t(z;h,n M ) = — , (2) 

Jo x^\l M x ■ i + {l-il M ) 

where h = H /100 km s" 1 Mpc -1 , and from now on the subscript denotes present day 
quantities. Note that for Qm = 1 the above expression reduces to the well known result 
for Einstein-de Sitter model (CDM, fi M = 1) for which t(z) = fif^l + z)~ 3 / 2 . As one 
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may conclude from the above equation, limits on the cosmological parameters Qm and H 
(or equivalently h), can be derived by fixing t(z) from observations. Note also that the age 
parameter, T = H t z , depends only on the product of the two quantities H and t z , which 
are usually estimated from completely independent methods (Alcaniz & Lima 1999). 

McCarthy et al. sample (2004) is constituted by 20 red galaxies on the interval 1.3 < z < 
2.2 selected from the Gemini Deep Deep Survey (GDDS). This sample consists of passively 
evolving galaxies whose most likely star formation evolution is that of a single burst of 
duration less than 0.1 Gyr (also consistent with Gyr), thereby meaning that the galaxies 
evolved passively since their first star burst formation. Another data set of interest is the 
Roche et al (2006) sample, which consists of 16 Extremely Red Galaxies (ERGs), formed 
by ellipticals, mergers and red spirals at 0.87 < z < 2.02. The authors have determined a 
deep spectroscopic survey of these galaxies on the Great Observatories Origins Deep Survey 
(GOODS)-South Field by using the Gemini multi-object spectrograph on the 8-m Gemini 
South Telescope. They also performed an age-dating analysis by fitting the spectra and 
nine-band photometry (BVIZJHK, plus 3.6/4.5-//m fluxes from Spitzer) of the ERGS with 
two-component models, consisting of passively evolving, old stellar populations combined 
with a younger continuously star forming component. The best-fitting mean ages for the 
old stellar populations range from 0.6 to 4.5 Gyr with an average value of 2.1 Gyr. The 
third set of data worked here is the Longhetti et al. (2007) sample, composed of nine field 
galaxies spectroscopically classified as early-types at 1.2 < z < 1.7. They result from a near- 
IR spectroscopic follow-up of a complete sample of bright ERGs selected from the Munich 
Near-IR Cluster Survey (MUNICS) that provides optical (B, V, R, I) and near-IR (J and K') 
photometry. The ages of the galaxies were estimated by making use of stellar populations 
with the degeneracy between the age of the best-fitting models and the star formation time 
scale being broken through a mass-weighting of the masses. As one may check, this sample 
provides the most accurate ages and the most restricting galaxies. 

In Fig. 1 we show a plot of the age as a function of the redshift for all objects present 
on the above quoted samples. However, since the galaxy formation is a random process it 
is not interesting to consider only the ages of all objects in order to constrain the model 
parameters. One basic reason comes from the fact that whether a universe model is capable 
of explaining the existence of a very old object for a given value of z, it explains more easily 
the age of the youngest objects at that redshift. With basis on this kind of criterion we have 
selected a more restrictive subsample involving only the oldest objects for a fixed redshift 
interval. The complete subsample chosen from the 3 set of data is constituted by 13 galaxies. 
As can be seen in the above quoted papers, some galaxies have asymmetric age uncertainties, 
so we have used the D'Agostini (2004) method to symmetrize it in order to perform a \ 2 
statistical analysis. 
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In Figs. 2 and 3 we show the selected subsample and the predictions of several ACDM 
models, in order to verify visually the concordance with the data and the effects of the 
different parameters (h and flu)- 

3. Statistical Analysis and Results 

Now, let us perform a x 2 fit over the h — Vt M plane. To begin with, we first consider 
that the actual ages of the objects at a given redshift are t z = t° bs + t inC) where t inc is the 
incubation time of the galaxy. This quantity is an estimate of the amount of time interval 
from the beginning of structure formation process in the Universe until the formation time 
(tf) of the object itself. It is also assumed that t inc varies slowly with the galaxy and redshift 
in our sample, and, in order to account for our ignorance about this kind of "nuisance" 
parameter, we associate a reasonable uncertainty, say, a tinc (Fowler 1986, Sandage 1993). In 
what follows, we consider that t inc = 0.8 ± 0.4 Gyr. Following standard lines, the maximum 
likelihood, C oc exp [— x 2 age ( z '-> P>^nc)/2], is determined by a x 2 statistics 

2 / I \ \ ^ P) tobs,i tinc) /„n 

Xage\ Z W) ~ ~2 ~7~Z2 ' \ 6 > 

l L obs,i Hnc 

where p = (h, Q m ) is the complete set of parameters, t b s ,i is the observational age for a 
specific galaxy, a tobs . is the uncertainty in the individual age, a tinc is the above assumed 
uncertainty of the incubation time. 

It now proves convenient both from a qualitative and a methodological viewpoint to 
consider first the old galaxy age test separately. Further, we present a joint analysis including 
the BAO signature from the SDSS catalog. It is worth noticing that a specific flat cosmology 
has not been fixed by hand in the analyzes below. 

3.1. Limits from Old Galaxies 

Let us now consider the selected 13 galaxies which constitute the more stringent old 
galaxies present in the original samples. For the sake of definiteness this is the sample 
adopted for deriving all the limits in this work. In principle, it is reasonable to think that 
old galaxies at high-z alone can not put very tight constraints on the pair (h,Q M ). One of 
the main reasons is that the galaxy formation is a random event, and, as such, it is hard to 
know the correct values of the incubation time for each galaxy. As it appears, the incubation 
time for some galaxies may be larger or smaller than the ones assumed here. 
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Fig. 1. — Plane Age- Redshift and total sample of galaxies. Green points correspond to the 
McCarthy et al. (2004) sample while black and red points represent the Roche et al. (2006) 
and Longhetti et al. (2007) samples, respectively. 
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Fig. 2. — Age-Redshift diagram (the effect of h). The data points correspond to the 13 
galaxies of our selected subsample (see the main text). The green and red filled circles stand 
for the observed age and the age plus the median of the incubation time (ti nc = 0.8Gyr) 
for each object. The dotted curves are the predictions of the cosmic concordance model 
(fl M = 0.3, f^A = 0.7) for some different values of the h parameter. As expected, for small 
values of h it is easier for the models explain the existence of old objects. 
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Fig. 3. — Age- Redshift diagram (the effect of Qm)- As before, the data points correspond 
to the 13 galaxies of the selected subsample with green and red filled circles standing for the 
observed age and the age plus the median of the incubation time {t inc = 0.8Gyr). Different 
from the previous figure, now we have plotted models with h = 0.72, as favored by the HST 
Key Project, for selected values of the density parameter. For smaller values of Qm the ages 
increase thereby accommodating the oldest selected objects. 
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In Fig. H] we show the contours of constant likelihood (68.3%, 95.4% and 99.7% c.l.) 
in the space parameter h — Qm for the age data alone. Note that all the parameter space 
for h and Qm is allowed, with a weak constraint Qm ^ 0.02 and h > 0.34. In this case, 
the best fit parameters derived are h = 1.11 and Qm = 0.09 with a Xmin = H-2, however, 
these values are meaningless because many others pairs of (h, CIm) are also allowed with 
considerable statistical meaning. For example, by fixing Qm = 0.3 we have h = 0.68 with 
Xmin = 11-4, and if VL M = 1.0 we have h = 0.67, and both cases are permitted with high 
degree of confidence. Clearly, one may conclude that an additional cosmological test (fixing 
JIm) is necessary in order to break the degeneracy on the (h, Qm) plane. 



3.2. Old Galaxies and BAO: A Joint Analysis 

As suggested earlier, more stringent constraints on the space parameter (Qm, h) can be 
derived by combining the selected sample of old galaxies with the BAO signature (Eisenstein 
et al. 2005). According to cold dark matter (CDM) picture of structure formation, the large- 
scale fluctuations have grown since z ~ 1000 by gravitational instability. The cosmological 
perturbations excited sound waves in the primeval plasma, thereby producing the acoustic 
peaks in the early universe. The peak detected (from a sample of 46748 luminous red 
galaxies selected from the SDSS Main Sample) is predicted to arise precisely at the measured 
scale of 100 h~ l Mpc. Basically, it is a consequence of the baryon acoustic oscillations in 
the primordial baryon-photon plasma prior to recombination. Let us now consider it as 
an additional cosmological test over the age redshift test. As it is widely known, such a 
measurement can be characterized by the dimensionless parameter 



Q 1/2 

A~ 



H{z if 



1/3 



— r(^* 



2/3 

= 0.469 ± 0.017, (4) 



where z* = 0.35 is the redshift at which the acoustic scale has been measured, and r(^) 
is the dimensionless comoving distance to z*. The quantity TC(z*) corresponds to H(z)/H , 
which, in a ACDM flat cosmology reads: 

n = [n M {i + zf + i-n M \ 1/2 (5) 



Note that the quantity given by (jl]) is independent of the Hubble constant and, as such, 
the BAO signature alone constrains only the Qm parameter. However, as can be seen from 
Fig. 3, the value of h is heavily dependent on Qm thereby leading to a tighter constraint 
in the distance scale. This property is very characteristic of the BAO signature thereby 
differentiating it from many others classical cosmological tests, like the gas mass fraction 




Fig. 4. — Confidence regions (68.3%, 95.4% and 99.7%c.l.) in the (fl m , h) plane provided 
by the Ages of 13 Old Galaxies at High Redshift. The best fit values are h = 1.11 and 
Q M = 0.09. 
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Fig. 5. — Contours in the VLm — h plane using the Age Redshift and BAO joint analysis. 
The contours correspond to 68.3%, 95.4% and 99.7% confidence levels. The best-fit model 
converges to h — 0.70 and f2 M = 0.28. 
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Fig. 6. — Likelihood function for the h parameter in a flat ACDM universe, from the age 
of old galaxies at high redshifts in combination with BAO signature. The shadow lines are 
cuts in the regions of 68.3% and 95.4% probability. We see that the region permitted is well 
constrained and in concordance with others studies (Freedman et al. 2001; Spergel et al. 
2006). 
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(Lima et al. 2003; Allen et al. 2004), luminosity distance (Peebles & Ratra 2003), the age 
of the universe (Alcaniz et al. 2003; Lima 2004; Cunha & Santos 2004; Jesus 2006), angular 
diameter versus redshift (Lima & Alcaniz 2000; 2005), and the Sunyaev-Zeldovich effect 
(Cunha et al. 2007). 

In Fig. [5j we show the confidence regions for the age redshift relation and BAO joint 
analysis. By comparing with Fig. HJ one may see how the BAO signature breaks the degen- 
eracy in the (f^M, h) plane. As it appears, the BAO test presents a striking orthogonality 
centered at Qm = 0.27 ± 0.03 with respect to the age redshift data as determined from old 
Galaxies at High Redshift. We find h = 0.71 ± 0.04 and xLn = n -38 at 68.3% (c.l.) and 
95.4% (c.l.), respectively, for 1 free parameter. An important lesson here is that the com- 
bination of old galaxies with BAO provides one more interesting approach to constrain the 
Hubble constant. In Fig. [HI by projecting over fl M , we have plotted the likelihood function 
for the h parameter based on our joint analysis involving the BAO signature plus the OHRG 
data. The dotted lines are cuts in the regions of 68.3% probability and 95.4%. 

As can be seen from Table [TJ the present results are in line with some recent analyzes 
based on different cosmological observations, like the one provided by the HST Project 
(Freedman et al. 2001) and the WMAP team (Spergel et al. 2007). Note, however, that it is 
just marginally compatible with the determination advocated by Sandage and collaborators 
(2006), a result obtained with basis on Type la Supernovae calibrated with Cepheid variables 
in nearby galaxies that hosted them. Jimenez et al. (2003) also obtained an independent 
estimate of the Hubble constant by computing dz/dt at z ~ for a sample of galaxies 
from Sloan Digital Sky Survey (SDSS). By assuming a ACDM concordance model, they 
obtained H$ = 69 ± 12 km s _1 Mpc -1 . More recently, Cunha, Marassi &; Lima (2007) have 
found h = 0.74+ 

ao3 an d — 0.28 ± 0.04 at 68% c.l. through a joint analysis involving 
the Sunyaev-Zel'dovich effect (SZe), X-ray surface brightness of galaxy clusters and baryon 

Table 1. Limits to H 



Method Reference 



CMB Spergel et al. 2006 (WMAP) 0.73 ± 0.03 

Cepheid Variables Freedman et al. 2001 (HST Project) 0.72 ± 0.08 

Age Redshift Jimenez et al. 2003 (SDSS) 0.69 ± 0.12 

SNe Ia/Cepheid Sandage et al. 2006 0.62 ± 0.013(rand.)±0.05(syst. 

SZE+BAO Cunha, Marassi & Lima 2007 0.74t°;^ 

Old Galaxies + BAO this letter 0.71 ± 0.04 



-14- 



acoustic oscillations. They have also assumed a ACDM flat cosmology. The present result 
is in perfect agreement with their analysis with a high degree of statistical confidence (both 
analyzes are compatible even at 68% c.L). Perhaps more important, the result is heavily 
dependent on the combination of the samples. For instance, if we have taken only the Roche 
et al. sample (the less restrictive one), after the joint analysis with BAO, the best fit increases 
toh = 0.99 (h = 0.77 for the McCarthy et al. data). 

It is also worth noticing that the best-fit scenario derived here, = 0.27 ± 0.03 and 
h = 0.71±0.04, corresponds to an accelerating Universe with q = —0.58, a total evolutionary 
age of t Q ~ 9.6/i -1 Gyr, and a transition redshift (from deceleration to acceleration) z t — 0.72. 

4. Conclusions 

Since the original work of Hubble, the precise determination of the distance scale, H , 
has been a recurrent problem in the development of physical cosmology. In this letter 
we have derived new constraints to the Hubble constant based on a statistical analysis of 
Ages of Old Galaxies at High Redshifts. Our analysis revealed that with basis only on 
the present age-redshift data of old galaxies it is not possible to obtain tight constraints 
on H due to the dependence of t(z) with the Vt M parameter (see Fig. 4). However, the 
degeneracy on the Vt M — h plane has been broken trough a joint analysis involving the baryon 
acoustic oscillation signature from the SDSS catalog (see Fig. 5). The Hubble constant was 
constrained to be h = 0.71 ± 0.04 (la). These limits were derived assuming an incubation 
time t inc = 0.8 ± 0.4Gyr and a flat ACDM scenario. 

As we have seen, the baryon acoustic oscillations (BAO) provide an interesting tool for 
constraining directly the mass density parameter, Qm, and, indirectly, it also improves the 
Hubble constant limits acquired from other cosmological techniques (like the age-redshift 
relation). Our Hubble constant estimation through a joint analysis involving old high red- 
shift galaxies plus BAO is in good agreement with independent estimates based on different 
cosmological observations, like the third year of the WMAP and the HST Key Project. 
Indirectly, such an agreement suggests that the incubation time we have assumed is quite 
realistic. It also reinforces the interest to the observational search of old galaxies at high 
redshifts in the near future, when more and larger samples, smaller statistical and systematic 
uncertainties can improve the limits on the present value of the Hubble parameter. 
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